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Abstract. The mean- field limit for the dynamics of bosons with random 
interactions is rigorously studied. It is shown that, for interactions that are 
almost surely bounded, the many-body quantum evolution can be replaced in 
the mean-field limit by a single particle nonlinear evolution that is described by 
the Hartree equation. This is an Egorov-type theorem for many-body quantum 
systems with random interactions. 



1. Introduction 

This work is a modest contribution to the mathematical theory of the mean- 
field limit for bosons with random interactions. There has been substantial de- 
velopments in the study of the mean-field dynamics of bosons with deterministic 
interactions. Early results were proven by Hepp in PQ, see also [2]. A differ- 
ent approach based on the reduced density matrix was developed in [3] and was 
substantially extended to more general potentials and to the derivation of the 
Gross-Pitaevskii equation in [I], [5], [6], [7J, [8]. Recently, a new approach was 
developed in [9], which gives convergence estimates in the mean-field limit that 
are uniform in Planck's constant h, see also |lUj . 

While the mean-field dynamics for bosons with deterministic interactions has 
attracted considerable interest, the question of the mean-field dynamics of bosons 
with random interactions has not been addressed, yet. Many-body bosonic sys- 
tems with random interactions are relevant to concrete physical systems, such 
as inhomogeneous nonlinear optical media, or Bose-Einstein experiments where 
irregular fluctuations in currents inside conductors close to the condensate induce 
via Feshbach resonances inhomogeneous interactions between the bosons, see [11] 
for a description of the latter; also [12] and references therein. Here, we give a 
simple recepy for extending the deterministic mean-field analysis to the case of 
random interactions (and in the presence of a random potential). 

1.1. The model. Consider the probability triple (H,JF, P), such that the prob- 
ability space Q has a generic point uj and is endowed with measure /i. Define on 
this space the random field 

v(x,u) : M 3 xfi^R, 
l 
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such that v is measurable in x £ R 3 and u G f2, and is almost surely in L°°(IR 3 ), 
i.e. there exists f2 C f2 such that /x(fio) = 1 an d, for all u G f2 , u (•, u) G L°°(IR 3 ). 
A concrete example of v that satisfies the above conditions is v(x,u) = v\(x) + 
V2(x, u), such that v\ G L°° and ^2 is Gaussian with finite mean and variance. 
For a measurable and integrable function / on Q, we define the expectation value 
of / as 

E(/) := J f{u)n{dw). 
We consider the iV-body random Schrodinger operator 

N 1 

(1) H* = H? := + ± Yl "0* 

i=l l<i<j<N 

where A = z2j=i ^ * s the 3-dimensional Laplacian and u G Q. Here, we work 

in units where Planck's constant H — 1 and the mass of each particle is m = \. 
We note that the analysis below is uniform in h. The Hamiltonian H N acts on 
the Hilbert space 7i N := L S (M. 3N ), the symmetrization of L 2 (M 37V ), which is the 
space of pure states for a system of N nonrelativistic bosons. 

The quantum dynamics of the iV-body system is described by the Schrodinger 
equation 

(2) idtV N {t)=H N * N (t), 

with an initial condition ty N (t = 0) = * 7V '° G L|(IR 3Ar ). 

Together with the dynamics defined above, the iV-body system is described 
by a kinematical algebra of "observables" . For p < N, a p-particle observable is 
described by an operator G B(TC^ P '), where B(Jiy p ') is the algebra of bounded 
operators on Tv?) — L S (B? P ). By the nuclear theorem, one can associate with 
a® a tempered distribution kernel in iS'(IR 3p x R 3p ), a^(x±, . . . , x v \ y±, . . . , y p ) : — 
a ip) (X p ;Y p ), such that 

(3) (a^V W )PQ= f a <*\X p ;Y p )tpV)(Y p )dY p 

where ip^(Y p ) G L|(IR 3p ). We associate to a^' an operator A N (a^) acting on 
that is given by 

(4) (A N (a^)( Xl ,---,x N ) = ™_ (P ga W ® Ii»-P)p 8 *)(xi, ...,x N ), 

where ^(xi, . . . , Xjv) G L|(IR 3Ar ) and Pg is the projection onto the symmetric 
subspace L 2 S (R 3N ) of L 2 (R 3N ). It follows from (J3J) and flU) that the map 

A N : B{H {p) ) ^ B(H {N) ), l<p<N, 
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is linear, such that 

(5) II^(« (P) )II B (^))<II« (P) II B ( W (,)), 

A N (a {p) Y = A N (a {p) *). 

In the Heisenberg picture, the evolution of A^ G B(Ti.^) is given by 

(6) of (iW) := e iHNt A^e~ iHNt , t G R. 

Since w is almost surely bounded, i7 is almost surely self-adjoint on the sym- 
metrized Sobolev space Hg(M. 3N ), and hence the propagator e~ lH t G M, 
is almost surely unitary. Moreover, it follows from the fact that the pointwise 
limit of measurable functions is itself measurable, [13], and the Trotter product 
formula, [H], that 

is u;- measurable. 

We now introduce the classical evolution. The Hartree equation is given by 

(7) i^t = -A^ + («*|^t| 2 )V't > 

with the initial condition ipt=o — <ft G L 2 (R 3 ). It follows from Duhamel's formula 
for ^ and the fact that t> G L°° almost surely, that global solutions of (J7|) in L 2 
exist almost surely, such that H^tlU 2 = ||0||l 2 with probability 1, for all i 6 1, 
(see for example [TH] for the case when v G L°°). It also follows from Duhamel's 
formula that the random variable 

(®f =1 ^, A® ® p i=1 ^), A® g 

is uj- measurable. 

1.2. Statement of the main result. We are in a position to state the main 
result. 

Theorem 1. Given a^ p \ A N (a^) and af as above, suppose that the initial state 
of the N-body system is a normalized coherent (product) state ^ N,0 (xi, • • • , xn) = 
®? =1 (j){xi), <f> G L 2 (R 3 ). Then, for fixed t > 0, 

lim E((V N >°,a?(A N (aV>))y N > )) = E«®? =1 ^, a 00 ®f =1 

N— >oo 

where ipt satisfies the Hartree equation ffity with initial condition ipt=o = 4>- 

We note that the analysis below can be easily extended to study the mean-field 
dynamics of bosons in a random external potential that is almost surely smooth, 
polynomially bounded and positive, and to investigate the semi-classical limit of 
the dynamics under additional assumptions on the decay of the interaction, as 
in [9]. Furthermore, the analysis below can be applied "in toto" to extend the 
results of [10] and [16] to the case of random interactions. 



2. Proof of Theorem Q] 

The proof of Theorem [1] follows effectively from an application of the dominated 
convergence theorem, see [T3|, and Theorem 1.1 in [9]. In what follows, we drop 
the explicit dependence on the time t in the notation, since we fix it. 

Proof. We introduce the random variables 

x (p) ._ (pNfl ia N(A N (aM))V N ' ) 

and 

X (P) ■= (®Li^,a (p) ®U 
The claim of the theorem is equivalent to the statement 

(8) lim E(X^) = E(X (p) ). 

N— >oo 

We divide the proof of (jHJ) into several steps. 
Step 1. Uniform integrability. We want to show that 

(9) BmE(|X«|l |Jf „,^)=0 J 

uniformly in N G N. 

We have from (0) and the fact that the quantum time-evolution is almost surely 
unitary, that 

(10) |X^ P) | < \\a^\\ B{7iip)) <2\\a^\\ b(«Cp)) < °° ; almost surely, 
uniformly in N G N. For > 0, it follows from fflQj) that 

(11) 1^11 | X (p)| >/3 < l-^w I < 2||a^ ||b(^(p)) < oo, almost surely, 
uniformly in N G N. The dominated convergence theorem together with (fTTT) give 

S'iep 2. Mean-field limit with probability 1. It follows from the fact that the 
particle interaction v G L°° almost surely and Theorem 1.1 in [5] that, for fixed 
t > 0, 

(12) X$ N =?°XV) almost surely. 
Step 3. It follows from Fatou's lemma, [T3], and fjTD]) . that 

(13) E(|XW|) < liminfE(|X^|) < limsup E(|X^|) < 2||a< p > || B(w(p)) < oo, 

uniformly in N G N. We also have that 

\x^\i lx(P)l ^<\x^\, 

which together with (TlBl and the dominated convergence theorem, imply that 

(14) limE(|XW|l ^=0. 



Step 4- Convergence as N — > oo. We introduce the random variable 

y(p) ._ I v-(p) _ v- (p) l 

Note that it suffices to show that E(Y^) — > as iV — ► oo, from which (JHJ) follows 
by the triangular inequality. 

It follows from ([12]), Step 2, that 

(15) ^^4°° o almost surely. 
We decompose into two parts, 

V (p) _ V b),<P , v(p).>/3 
1 N — 1 N ' N ' 

where if ),</3 := ^l^w,^ and Y*^ := V^l,^,,^, for P > °" 

Since Y^' </3 < /3, (1151) together with the dominated convergence theorem 
imply that 

(16) lim E(y$ ),</3 ) = 0. 

N— »oo 

Furthermore, since 

r N s '\±\x(p)\>f3/2 + z \^n l 1 !^!^^' 
it follows from <Q and flH|) that 

(17) lim E(Y^- P ) = 0, 

uniformly in N G N. 

Given e > 0, ( 1TT1) implies that there exists a finite /5o > such that 

supE(y^ ) '- /3() ) < e/2. 

N 

Moreover, (1T61) implies that there exists a positive integer Nq such that, for all 
N>N Q , 

E (yb),</3o ) < e/2 _ 

It follows that 

E(4 P) ) = E(4 p) ' </3 °) + E(Y^°) < e 

for N > N Q . Therefore, E(Y^ p) ) N =?° 0. 
By the triangular inequality, 

|E(X&0) _ E(X^)| < E(F i {f ) ) ^4.°° 0, 

which gives the claim of the theorem. □ 
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